求解一维初值问题:特征线法与D'Alambert公式

#特征线法 #D_Alambert公式 #半无界问题 #对称开拓法 #奇延拓

我们求解维数 n=1 时的波方程初值问题

WaveCauchy(f,φ,ψ):={Lu=utta2x2u=f(x,t),(x,t)R×R+,u(x,0)=φ(x),xR,ut(x,0)=ψ(x),xR.

初值问题的求解

特征线法

注意到此时 =x,此时可作分解

L=t2a2Δ=(tax)(t+ax)

v:=(tax)uWaveCauchy(f,φ,ψ) 可化为如下的两个初值问题:

{utaux=v,u(x,0)=φ(x);{vt+avx=f,v(x,0)=ψ(x).

然后分别用一阶线性方程与特征线法 中提到的特征线法分别求解下式、上式即可.

D'Alembert 公式

上面的特征线法求解的结果为

u(x,t)=φ(x+at)+φ(xat)2+12axatx+atψ(ξ)dξ+12a0tdτxa(tτ)x+a(tτ)f(ξ,τ)dτ

f0,即原方程为齐次方程时,方程的解为

u(x,t)=φ(x+at)+φ(xat)2+12axatx+atψ(ξ)dξ

上式称为 D'Alembert 公式

D'Alembert 公式可表示为

()u(x,t)=F(x+at)+G(xat)

例如

F(s)=12φ(s)+12a0sψ(ξ)dξ,G(s)=12φ(s)12a0sψ(ξ)dξ

物理意义为左行波与右行波的叠加.

另一方面,注意到形如 () 的解均为一维波方程(不考虑初值条件)utta2uxx=0 的解,而 () 中出现的两项分别为 utaux=0ut+aux=0 的解.

相容性条件:形式解为古典解的条件

φC2 的,ψC1 的,fR×[0,+)C1 的.

此时解 uC2 的,从而是古典解.

推论:解的奇、偶、周期保持

如果 f,φ,ψ 都是奇/偶/周期为 l 的函数,则解也为奇/偶/周期为 l 的函数.

半无界问题的求解

方程组

{Lu=f(x,t),in R+×(0,+)u=φ, ut=ψ,on R+×{t=0}u=g,on{x=0}×(0,+)

称为半无界问题.求解半无界问题的想法是将定解问题的初值和非齐次项延拓到整个实数轴上,从而转化为 Cauchy 初值问题求解.

齐次边值半无界问题:对称开拓法

g0 时,称半无界问题为齐次半无界问题,即

{Lu=f(x,t),in R+×(0,+)u=φ, ut=ψ,on R+×{t=0}u=0,on{x=0}×(0,+)

考虑对 x 作奇延拓:

φφ¯,ψψ¯,ff¯

这里 ¯(x,t)={(x,t),x0(x,t),x<0 .则原题化为 Cauchy 初值问题

WaveCauchy(f¯,φ¯,ψ¯).

方程的解 u¯ 依推论为奇函数.由 Cauchy 初值问题解的表达式知

u¯(x,t)=φ¯(x+at)+φ¯(xat)2+12xatx+atψ¯(ξ)dξ+12a0tdτxa(tτ)x+a(tτ)f¯(ξ,τ)dξ.

进一步化为:

xat 时,

u(x,t)=φ(x+at)+φ(xat)2+12axatx+atψ(ξ)dξ+12a0tdτxa(tτ)x+a(tτ)f(ξ,τ)dξ.

x<at 时,

u(x,t)=φ(x+at)+φ(atx)2+12aatxx+atψ(ξ)dξ+12a(txatdτxa(tτ)x+a(tτ)f(ξ,τ)dξ+0txadτa(tτ)xx+a(tτ)f(ξ,τ)dξ).

相容性条件

若各项满足下列条件:

  1. φC2φ(0)=0
  2. ψC1ψ(0)=0
  3. fC1a2ψ(0)+f(0,0)=0.

则上面给出的 uC2 ,为齐次半无界问题的解.

说明

u(0,0) 连续推 φ(0)=0ut(0,0) 连续推 φ(0)=0uC2 且在 (0,0) 处适合方程推 a2ψ(0)+f(0,0)=0. 注意应用边界 {x=0}×(0,+)u=0 的条件.

非齐次边值半无界问题

一般地,对非齐次半无界问题

{Lu=f(x,t),in R+×(0,+)u=φ, ut=ψ,on R+×{t=0}u=g,on{x=0}×(0,+)

作变量替换:

u(x,t)=v(x,t)+g(t)

代入得

{Lv=f(x,t)g(t),in R+×(0,+)v=φg(0), vt=ψg(0),on R+×{t=0}v=0,on{x=0}×(0,+)

从而,该非齐次边值半无界问题等价于新的齐次边值半无界问题.

第二边值条件非齐次边值问题

若改为第二边值条件,变为:

{Lu=f(x,t),in R+×(0,+)u=φ, ut=ψ,on R+×{t=0}ux=g,on {x=0}×(0,+)

则可作函数变换

u(x,t)=xg(t)+v(x,t)

将边值问题转化为齐次边值条件

vx=0,on {x=0}×(0,+).

再由对称开拓法(偶延拓)即可.求解多维初值问 M